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Abstract 

Isoscalar Cooper pairing in isospin asymmetric nuclear matter occurs between 
states populating two distinct Fermi surfaces, each for neutrons and protons. 
The transition from a BCS-like to the normal (unpaired) state, as the isospin 
asymmetry is increased, is intervened by superconducting phases which spon- 
taneously break translational and rotational symmetries. One possibility is 
the formation of a condensate with a periodic crystallinelike structure where 
Cooper pairs carry net momentum (the nuclear Larkin-Ovchinnikov-Fulde- 
Ferrell-phase) . Alternatively, perturbations of the Fermi surfaces away from 
spherical symmetry allow for minima in the condensate free energy which cor- 
respond to a states with quadrupole deformations of Fermi surfaces and zero 
momentum of the Cooper pairs. In a combined treatment of these phases we 
show that, although the Cooper pairing with finite momentum might arise as 
a local minimum, the lowest energy state features are deformed Fermi surfaces 
and Cooper pairs with vanishing total momentum. 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The pairing properties of nuclear systems - finite nuclei and the bulk nuclear matter - 
play an important role in the physical manifestations of these systems. In most cases of 
interest the pairing occurs at finite isospin asymmetry. For example, charge neutral bulk 
nuclear matter with an admixture of Icptons in 3 equilibrium implies proton abundancies at 
about 40% in supernova matter and 10% in neutron star matter. The dominant partial-wave 
channel for the pairing depends on the density, temperature, and the isospin asymmetry in 
general. For large enough asymmetries the matter is paired in the isospin triplet channels 
- the ^Sq channel at low densities [1-11] and ^P2-^-^2 or channels at high densities 
[12-20]. For weakly asymmetric systems the isospin sing let attractive ^Si-^Di [22-27] and 
channels [28,29] dominate the pairing interaction at low and high densities, respectively. 
In bulk nuclear matter of compact stars the asymmetries are most likely too large to support 
isospin singlet pairing. However the dilute nuclear matter in supernovas and the low-density 
tails of the exotic nuclei can support ^5'i-^£)i-channel pairing, since in general the isospin 
asymmetry is less effective in destroying the pair correlations in dilute matter [25]. 

The mechanism of the suppression of the BCS pairing of the kind found in asymmetric 
nuclear matter was initially studied in metallic superconductors in a spin-polarizing magnetic 
field [30-33] . Since a homogenous magnetic field is screened beyond the London penetration 
depth in the bulk of a superconductor, in real materials the electron spin polarization is 
induced by low-density paramagnetic impurities, which interact with an electron by spin 
nonconserving forces. The effect of collision induced polarization is commonly modeled by 
an average polarizing field which implies an asymmetry in the populations of spin-up and 
-down quasiparticles. The resulting electron Pauli paramagnetism suppresses the 5"- wave 
pairing when the Fermi surfaces of the electrons with spin-up and -down are apart on a 
scale of the order of the energy gap in the quasiparticle spectrum. A close analogy to the 
ordinary superconductors can be seen in the effect of the large B fields in the neutron stars 
on the S'-wave pairing of neutrons. Our general discussion below also applies to the pairing 
in spin-polarized neutron matter in the ^5*0 channel, although the numerical computations 
are carried out for the isoscalar pairing in asymmetric matter. 

The understanding of the isospin singlet pairing as a function of the isospin asymmetry 
requires an extension of the standard BCS theory to describe phases with spontaneously bro- 
ken space symmetries. Larkin and Ovchinnikov [34] and Fulde-Ferrell [35] (LOFF) argued, 
in the context of metallic superconductors, that a transition from the BCS paired to the 
normal state occurs via a phase where the Cooper pairs carry a finite total momentum. The 
isospin polarization in nuclear matter with pairing in the ^Si-^Di channel leads to a nuclear 
LOFF phase for a range of isospin asymmetries [36,37]. Unlike the metallic superconductors 
where the ratio of the pairing gap to the average chemical potential is of the order 10~^, 
the nuclear LOFF phase emerges in the strong coupling regime, where this ratio is of the 
order 0.1. The spatially periodic structure of the nuclear LOFF phase corresponding to 
oscillations of the condensate wave function are on the scales of the order of several fm. 
Note that, independent of the detailed structure of the lattice, the LOFF phase breaks both 
the translational and rotational invariance of the original BCS state. The formation of the 
LOFF phase is a robust feature of the fermionic systems with broken time-reversal symme- 
try and has been studied in detail in the fiavor asymmetric high-density QCD (crystalline 
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color superconductivity) [38-45]. 

The evolution of the condensate from the BCS to paired states with broken space sym- 
metries, as the isospin asymmetry is increased, can take a different path if we give up the 
assumption that the intrinsically homogeneous matter features spherical Fermi surfaces [46] . 
As shown in ref. [46], for quadrupole deformed Fermi surfaces (DFS) the BCS condensate 
evolves to a state with spontaneously broken rotational symmetry [in group theory terms 
the 0(3) symmetry breaks down to 0(2)]. Such a state is still axially symmetric and, unlike 
the LOFF phase, preserves the translational symmetry. We shall refer to the superconduct- 
ing state with deformed Fermi surfaces as the DFS phase with the understanding that the 
expansion describing the deformation of Fermi surfaces is truncated at quadrupole order 
(the higher order terms in the multiple expansion have not been studied so far). 

Why do the phases with broken space symmetries dominate the BCS state? The BCS 
quasiparticle spectrum (for homogeneous systems) is isotropic; when the polarizing field 
drives apart the Fermi surfaces of fermions with spin/isospin up and down the phase space 
overlap is lost, hence the pairing is suppressed. The finite momentum of the Cooper pairs 
of the LOFF phase implies an anisotropic quasiparticle spectrum; the magnitude of the 
anisotropy is controlled by the value of the net momentum of the pairs which is treated as 
a parameter. This new degree of freedom, which can be viewed as an additional variational 
parameter for the minimization of the ground state of the system, (partially) overcomes the 
phase-space loss caused by the departure of the Fermi surfaces from the perfectly overlapping 
configuration. In effect the gain in the pairing energy through increase of the pairing field 
dominates the increase of the kinetic energy. To motivate the DFS phase, note that formally 
the LOFF spectrum can be viewed as a dipole (oc Pi{x)) perturbation of the spherically 
symmetrical BCS spectrum, where Pi{x) are the Legendre polynomials, and x is the cosine 
of the angle between the particle momentum and the total momentum of the Cooper pair. 
The DFS phase assumes an expansion of the quasiparticle spectrum up to the next-to- 
leading order, quadrupole (oc P2{x)) deformation; it also assumed that the total momentum 
of the Cooper pairs is zero, i.e., there are quadrupole deformations only. Although the 
perturbations from the spherical symmetry cost kinetic energy, they are compensated by 
the gain in the potential energy. For a range of asymmetries the DFS phase becomes the 
ground state of the paired system [46] . 

The purpose of the present paper is a combined treatment of the LOFF and the DFS 
phases in bulk isospin asymmetric nuclear matter within the same model as that in Refs. 
[36,46]. The model assumes a pairing force given by a phase-shift equivalent (i.e., realistic) 
interaction and nonrenormalized quasiparticle spectrum. The effect of the mean field in the 
Brueckner-Hartree-Fock approximation is to reduce the range of isospin asymmetries where 
the pairing exists [24]. We shall follow the strategy of disentangling the effects of the mean 
field from the pairing since already at the level of approximations of Refs. [36,46] the LOFF 
and DFS phases show a complicated behavior. The renormalization of the quasiparticle 
spectrum for the pairing interaction is a problem of its own which has not been solved to 
date. Issues such as the consistency between the vertex and propagator renormalization and 
fulfillment of the spectral sum rules, etc. need further attention. 

In Section 2 we derive the BCS equations, which include the effects of the finite mo- 
mentum of the Cooper pairs and Fermi surface deformations, within the finite-temperature 
real-time Green's functions formalism. The numerical solutions of these equations are pre- 
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sented in Section 3, where we discTiss the phase diagram of the combined LOFF and DFS 
phases and their properties at finite temperatures. Section 4 contains a summary of the 
resuhs and an outlook. 



II. FORMALISM 



The starting point of the model is the Hamiltonian of fermions interacting via two-body 
forces: 



H 



d'^rVip^f, t, a)Vip{f, t, a) 



a/3 



where ip{f,t,a) are second quantized operators in the Heisenberg representation, Va,fi{r — 
r')6{t — t') is the space and time local bare interaction, and a, (5 stand for discrete quantum 
numbers (spin, isospin, etc.). We choose to formulate the finite-temperature pairing theory 
in terms of nonequilibrium real-time Green's functions. The single-particle propagator in 
this formalism is combined in a 2 x 2 matrix. 



G{Xi,X2) 



G''{xi,X2) G^{xi,X2) 
G>{XI,X2) G''(xi,X2) 



-i{Ti){xi)^\x2)) i{'4^\x2) '4^{xi)) 

-i{^{x,)li^\x2)) -i{f'lP{x,)'lP^{x2)) 



(2) 



where x = {f,t,a), and the averaging is over an arbitrary nonequilibrium state of the 
system; T and T are chronological and antichronological time ordering operators for the 
Nambu spinors 



— iIj^ (^x)lp{x) -0 



ij^x) 



(3) 



The Nambu spinors span the particle-hole space and satisfy the fermionic equal-time anti- 
commutation relations 



{M^,t),'iPl{x',t)} = 5^f,S'{x-x'), {Mx,t),Mx',t)} = 0, {M^,t)\i/;lix',t)} = 0. 

Starting from the equations of motion for the Heisenberg operators and the Hamiltonian (1) 
the Martin- Schwinger hierarchy can be constructed for a matrix Green's function involving 
an arbitrary number of field operators. For our purposes in is sufficient to truncate the 
hierarchy at the n = 1 level. The equation of motion of a one-particle Greens function is 
then given by the time-dependent Dyson equation 



t J,t/ 



G(o)(a:i) 0G{xi,X2) = (TzS{xi - X2) -V{xi,X6;X5,X4) ^ G{X5,X4,X2,X. 
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(4) 



where the subscript (0) refers to the free propagator, (g) stands for summation/integration 
over the repeated discrete/continuous variables, CTj are the components of the vector of the 
Pauli matrices, and the superscript -|- means an infinitesimal increment in the time argument. 
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The hierarchy at a given level can be decoupled (formally) by introducing the self-energy 
matrix: 

G(o](xi) (g) G{xi,X2) = azS{xi - X2) + X3) G{x3, X2), (5) 
which is defined as 

^{xi,X2) = —V_{xi,Xq;X5,X4) ^ G{x5,X4,Xi,X'^)G {Xt,X2). (6) 

The first of these equations defines the quasiparticle spectrum of the superconducting state, 
the second, the self-energies, e.g., the gap function. To close the set of exact equations we 
need to specify the approximation to the two-particle Green's function. The BCS theory fol- 
lows in the Hartree approximation, where the two-particle Green's function is approximated 
by a product of single-particle Green's functions. Note that the anomalous contributions of 
the type {ip{xi)ip{x2)) and ('(/'^(xi)-(/'^(x2)) are automatically built in in the Hartree approx- 
imation. This can be verified from the explicit form of the propagator in Eq. (2), when the 
spinors (3) are substituted. Each element of the 2x2 matrix (2) is in turn a 2 x 2 matrix 
in the particle-hole space: 

r<j _( <^^(^i'^2) F^{xi,X2)\ 

^''-\F'^^{X,,X2) G\X,,X2))' 

where j = c,a,>,<. The matrix form of the self-energies in the particle-hole space is 
analogous to that of the propagators: 

\T,>{xi,X2) T,''{xi,X2) J \ATJ(a;i,X2) E {xi,X2) J 

where the off-diagonal elements of the 'E^{xi,X2) matrix correspond to the pairing ampli- 
tude. Note that the inverse free-particle propagator is diagonal in the particle-hole, j and 
spin/isospin spaces: 

<^(0)"/3 = + VV2"^a - /^a) (9) 

The Dyson Eq. (5) is a 16 x 16 matrix equation in general if we include the spin and isospin. 
The size of the matrix to be diagonalized can be reduced in several successive steps. Here 
we note that in the time-independent stationary limit the number of the propagators can 
be reduced by applying the unitary transformation U — {1 + iay)/y/2 to Eqs. (5) and (6), 
in the j space, and introducing the retarded/advanced propagators by the relations 

^^'"^(^1, X2) = G%Xi, X2) - G<'>{xi, X2) = G'>'<(xi, X2) - ^"(xi, X2). (10) 

Further reductions will be carried out in the next section. The Dyson equation, e.g., for the 
retarded propagator is invariant under the rotation in the j space affected by the transfor- 
mation above, i.e., 

G^{xi, X2) = G'fo)(xi, X2) + Gfo)(xi, X2) ® S^(xi, X3) ® G^{x3, X2); (11) 

the retarded self-energies are related to the matrix components in the j space by relations 
analogous to Eq. (10). The integral Eq. (11) for the retarded propagators can be solved 
without specifying the approximation to (and hence the form of) the self-energies, in the 
quasiclassical approximation. 
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A. Quasiparticle spectrum 



When the characteristic length scales of the spatial variations of the macroscopic con- 
densate are much larger than the inverse of the momenta involved in the problem ~ Pf, 
where is the Fermi momentum, the two-point correlation functions can be approximated 
by their quasiclassical counterparts. Going over to the center of mass X — {xi + X2)/2 and 
relative ^ = Xi — X2 coordinates in the two-point functions and carrying a Fourier trans- 
form with respect to the relative coordinates we arrive at the mixed representation for these 
functions, e.g., 

G<(p, X)^l e'P^ G<(x + ^,X- d\ (12) 

where p = (a;,p) are the relative frequency and momentum, and X = {R,T). Since the 
variations of the propagators and self-energies are slow on the scales of the order of R, 
keeping the leading order terms in the gradient expansion is accurate to order ~ OippR)- 
The quasiclassical counterpart of the Dyson Eq. (11) is 



where 



e% = [(P/2 ± pf /2m« - ± eJ^+^(^, P/2 ± p) - eJ(-)(^, P/2 ± p). 



and i is a unit matrix in the particle-hole space. Here we have defined the symmetric and 
antisymmetric combinations of the particle/hole retarded self-energies 



Eff )(u;,p) ^ ^ E^,(u;,p) + E" (u;,p) , (14) 



-'a/3 



^R 



^R (, , ^^ 



The dependence on the ccntcr-of-mass time is dropped in the above expressions in the 
stationary limit. The retarded self-energies now can be expanded in Legendre polynomials: 

oo 

^at^~\^,P) = Y.^a&l~\^.P)Pl{G0S 6), (16) 
1=0 

where 9 is the angle between the vectors p and P. In practice the expansion will be truncated 
at order 1 — 2; the I — term renormalizes the chemical potential, as discussed below. For 
translationally invariant interactions the odd terms do not contribute in the expansion (16); 
the interactions are assumed to be timelocal. For translationally noninvariant interactions 
the effect of the / = 1 term would be a renormalization of the dipole deformation in the 
spectrum due to the kinetic energy terms. It can be included in the renormalization of the 
particle mass entering the anisotropic contribution to the kinetic energy, although it would 
not contribute to the net isotropic increase of the kinetic energy due to the finite momentum. 
Since the total momentum is treated as a variational parameter in the problem, the effect 
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of the sclf-cnergy renormalization at / = 1 can be accounted for by a simple rescaling of this 
parameter; therefore we anticipate that the / = 1 terms would not change the results below. 

For spin and isospin conserving forces the normal Green's functions and self-energies are 
diagonal in the spin and isospin spaces. It is sufficient to consider the anomalous propagators, 
e.g., in the isospin space, since the resulting spin structure, for 5"- wave interactions, is 
uniquely determined for each isospin combination. If we restrict ourselves to the neutron- 
proton pairing in the "^Si-^Di channel, which is justified when A„„, App <^ A„p (subscripts 
n and p refer to protons and neutrons) then A^p = cr^-A. The spectrum, in this case, takes 
the form 



(17) 



where the symmetric and asymmetric parts of the spectrum (which are even and odd with 
respect to the time-reversal symmetry) are defined as 



Ea 



2(^+-^-)- 



(18) 



The eigenvalues (17) are valid for an arbitrary approximation to the self-energies. Further, 
we shall approximate the symmetric and antisymmetric parts of the spectra using (i) the 
quasiparticle and (ii) the momentum-independent self-energy shift approximations. The first 
approximation implies 



''•^"/^^(p,ep)^ReSX/+)(p,6,) 



d 



P. + a^^^^"/^ 



p,(^p) {p-pf); (19) 



Pf 



the second approximation keeps the leading order (constant) term in the expansion (19). 
Within these approximations we obtain 



Es — 1 — i/ie + SeS/j,) cos^ 9, 

8m 2m 

Pp 

Ea = — cos 9 — 6fj, + (fj,6e + e6fj.) cos^ 9, 
2m 

where we defined the average chemical potential and the conformal deformation 



1 r 1 1 



e = — Re 
4/i 



^-{Pf) + %{Pf) 
s5-)(p^)+Ey(p^)' 



(20) 
(21) 

(22) 
(23) 



as the / = and / = 2 contributions from the isospin symmetric part of the quasiparticle 
spectra; similarly the relative chemical potential and the relative deformation are defined as 
the / = and / = 2 contributions from the isospin antisymmetric part of the quasiparticle 
spectra: 



1 1 

ARe[E5+)(p^) 



6e 



>,2 



(pf) 



(24) 
(25) 
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where we used the abbreviation 



n±{PF) = Re 



n/p 



(26) 



The self-energies on the right-hand side of Eqs. (22)- (25) are functionals of the pair- 
ing field, i.e., they are nonzero even in the case where all the interactions apart from the 
pairing interaction are switched off. They depend on the quasiparticle momentum via the 
momentum dependence of the gap function; if the latter is approximated by its value at 
the Fermi momentum the self-energies become constants. The values of the average and 
relative chemical potentials and 5^ can be adjusted to reproduce the particle density at a 
given temperature and isospin asymmetry. The conformal and relative deformations e and 
6e along with the finite momentum of the pairs P are treated as variational parameters to 
be determined from the ground state of the superconducting system. 

The solution of the Dyson Eq. (13) for the propagators can be written in terms of the 
eigenvalues as 



G^/,{u;,p) = + (27) 

F''{u,p) = UpVp ( ^— ^-—-) , (28) 

\ a; — + 17] u — + irj I 



where the Bogolyubov amplitudes are 

ul^-^ , vl^- , 29 

2 2y^E| + |A|2 " 2 2y^^| + |A|2 

The remainder Green's functions can be reconstructed using the spectral representation of 
the retarded Green's functions 



G„/,(-,p) = y_c^a;-^^-^, (30) 
F>,p)=r..'^i^, (31) 

J-oo UJ — UO' + IT] 

where A„/p(c(j,p) and B{uj,p) are the spectral functions. The Kadanoff-Baym ansatz [47] 
provides the link between the the retarded components and the remainder Green's functions 
in the j space. We have 

G>;;(c^,p)=A/^(a;,p)/>'<(u;), (32) 
F>'<{u;,p)^B{u;,p)f>'<{u;), (33) 

where, in equilibrium, the Wigner distribution functions of the Kadanoff-Baym ansatz [47] 
reduce to the Fermi distribution function f^{co) — [exp{/3cu) + and f^{uj) = [1 — /^(o;)]; 
here f3 is the inverse temperature. Note that the relations above are exact in the equilibrium 
limit. The causal and acausal Green's functions are then obtained from the relations (10). 
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B. The gap equation 

The BCS gap equation follows from the Hartree approximation for the two-particle 
Green's function in Eq. (4): 

G{xi, X2] X3, X4) = G{xi, X2) G{x3, X4), (34) 

which implies, according to Eq. (6), 

E(xi, X2) = -Y.{xi,X3; X2, X4) (g) G{x4, x^). (35) 

Upon applying the quasiclassical approximation to the above equation in the case of time 
local interaction, we obtain the retarded (off diagonal in the particle-hole space) self-energy 

A«(P,P) = / ^^V{p,p')F<{u;',p',n (36) 
or using the relation F^{u,p) = —2lmF{uj,p)f{u), which follows from Eqs. (31) and (33), 
A«(p,P) = -/ ^^V{p,p')lmF\u;',p',P)f<{u;'). (37) 

Further progress requires partial-wave decomposition of the interaction, which can be per- 
formed after angle averaging the remainder functions on the right-hand side of Eq. (37). 
The result of this procedure is 

Af(p,P) = - E / ^^KP,pO^|fi=#=([/<K) (38) 
{2ny 2^e| + A(p', P)2 

where (...) denotes the average over the angle between the relative and total momenta and 
A(]9, P)^ = Ao(p, P)^ -I- A2(]?, P)^ is the angle averaged gap. Here the pairing interaction 
is approximated by the bare neutron-proton interaction V{p,p') in the ^Si-^Di channel. 
As discussed above the average and relative chemical potentials can be fixed by adjusting 
the pairing field to reproduce the matter density p = Pn + Pp and the isospin asymmetry 
a — {pn — pp) I p. The total momentum, the relative and conformal deformations, are treated 
as variational parameters to be determined from the ground state energy of the system. 
The corresponding expression for partial densities of protons and neutrons is provided by 
the relation 

^^G'n/p(^>P)> (39) 

or, using Eqs. (30) and (32), 

Pnip = -2E/ ^ImG„/p(^,p)/(u;) = E/ {t.J/<(a;±) +^J/<(a;^)} , (40) 

where a stands for quasiparticle spin and the second relation follows in the quasiparticle 
approximation. 
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C. Thermodynamics 



The phase diagram of the paired state at fixed finite temperature and density can be 
obtained from the free energy in the mean field approximation. The free energy is given by 
the thermodynamic relation 

T\p,p = U~^-'S. (41) 

where U is the internal energy and S is the entropy. A thermodynamically stable paired 
state minimizes the difference of the free energies of the superconducting and normal states, 
= J^s — ^N- The mean field entropy of the superfiuid state is 

55--2^BE{r(^+)lnr(c^+) + r(^+)ln/>(^+) + /<(^-)ln/<(a;_) + />(^_)ln/>(^_)}, (42) 

where /c^ is the Boltzmann constant and the sum is over the momentum states in the 
quasiparticle approximation. The corresponding expression for the entropy of the normal 
state Sn is obtained by taking the limit A — > in Eq. (42). The internal energy of the 
superconducting state in the mean field approximation is 

^ ^ / + '~''^^^^] ^ ¥ ^''^^'^'^ ^i^p)Mp')\ (43) 

where the normal and superconducting occupation probabilities are defined as 



„/^(p) = ulf<{u^) + v'J<{u^), u{p) = u,v, [/<(a;+) - /<(a;_)] . (44) 

The first term in Eq. (43) includes the kinetic energy of quasiparticles which is a functional 
of the pairing gap. In the normal state it reduces to the kinetic energy of noninteract- 
ing quasiparticles. The second term describes the BCS mean-field interaction among the 
particles in the condensate and vanishes in the normal state. 



III. RESULTS 

The nuclear LOFF and DFS phases were studied numerically using the Paris nucleon- 
nucleon interaction. Our results did not depend on this choice of the interaction, since there 
are no significant deviations in the scattering phase shifts in the ^5*1 — ^ Di channel among 
various realistic potentials and these reproduce the properties of the deuteron at the same 
level of accuracy. The net density of the matter was fixed at the empirical saturation density 
of the nuclear matter ps = 0.17 fm^'^ At this density the typical energies of the nucleons 
are below the threshold laboratory energy E^ab — 220 MeV at which the interaction in the 
^Si —^Di channel becomes repulsive. Lower densities could provide a more realistic (from the 
physical point of view) setting, however, they have the disadvantage that the effects of the 
Bose-Einstein condensation of deuterons may start to play a role [25]. The computations 
were carried out for several fixed values of the density asymmetry a that are above the 
critical values at which the nuclear LOFF and DFS phases set in (see Refs. [36,46]). The 
qualitative behavior of quantities of interest is generic for different values of a, therefore, 
we shall show the results for a fixed value a — 0.35 at which the LOFF and DFS phases 
(separately) dominate the ordinary BCS state. 
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FIGURES 




FIG. 1. The pairing gap as a function of a Cooper pair momentum P in units of Fermi mo- 
mentum and the relative deformation parameter 6e. The density asymmetry is fixed at the value 
a = 0.35. The temperature and the density are T = 3 MeV and ps = 0.17 fm^^, respectively. The 
gap is normalized to its value Aqo in the asymmetric BCS state. 

The results of the combined treatment of the LOFF and DFS phases will be presented in two 
steps: in the next section we fix temperature and vary simultaneously the net momentum of 
the pairs P (LOFF phase) and the relative deformation Se (DFS phase). In the subsequent 
section we study the temperature dependence of these phases at either constant P or 5e. 

A. LOFF versus DFS 

Figure 1 shows the pairing gap as a function of the total momentum of Cooper pairs P and 
the relative deformation 6e for vanishing conformal deformation (e = 0). The temperature is 
fixed at the value T — 3 MeV. The ratios of the temperature over the Fermi energy and the 
pairing gap in the symmetric state are 0.1 and 0.4 respectively, i.e., they correspond to the 
low-temperature regime. The gap is normalized to its value Ago = A(P — 0,5e = 0) in the 
rotationally/translationally invariant state at a = 0.35. Although a = {pn — Pp)/{Pn + Pp) 
changes in the interval [—1; 1] in general, the symmetry of the equations with respect to 
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the indices labeling the species reduces the range of a to [0; 1]. The relative deformation 
obviously is not bounded and can assume both positive and negative values. Positive values 
of 5e imply an oblate deformation for the Fermi surface of neutrons and a prolate deformation 
for the Fermi surface of protons. The behavior of the pairing gap can be understood by 
examining the effect of the symmetric and antisymmetric parts of the quasiparticle spectrum 
on the gap Eq. (38). The antisymmetric part Ea appears only in the Fermi distribution 
functions. If Ea = one recovers the BCS limit, where the Fermi distribution functions 
become identical, i.e., the Fermi surfaces perfectly coincide. In this limit the phase-space 
overlap between the paired states is maximal. Consider the effect of a finite Ea when 
P — — Se. It increases/decreases the energies of the isospin up/down particles; the shift 
in the Fermi-levels of the isopin up/down particles samples different momentum regions 
of the phase-space in the kernel of the gap equation. This blocking effect or phase-space 
decoherence in the pair states in turn reduces the magnitude of the gap. In particular this 
allows for pairs with small Ea- Switching on a finite P and/or Se > acts to restore partially 
the phase space coherence of the isospin up/down states when 5/1 > max{Pp/2m, //5e} [see 
Eq. (21)]. Therefore, the magnitude of the gap increases in this case, i.e. the phase-space 
overlap of the Fermi surfaces is (partially) restored. Increasing P and/or 6e > further acts 
in a manner similar to pure energy shift Sjj, - the decoherence increases and the pairing is 
eventually destroyed. This picture is seen in Fig. 1 in the region Se > 0. Small perturbations 
in P and Se from the asymmetric BCS state {P — — Se but S/j, ^ 0) increase the pairing 
gap. For large perturbation in either P or Se > ^ the gap vanishes as the decoherence 
increases. The combined effect of the finite momentum and the quadrupole deformations is 
seen in the minimum plateau in the region of small P and 5e, followed by a maximum in the 
gap with increasing P and/or 5e, and a rapid falloff beyond the maximum. The behavior 
of the gap in the limits of pure LOFF and DFS phases is the same as the one found in the 
previous work [36,46]. 

In the case of 5e < the ^iSe term in Eq. (21) combines with the chemical potential shift 
S[i. To satisfy a given particle number asymmetry there are now two parameters (while there 
was only one, for positive Se). As a result with increasing deformation and hence \iSe (the 
changes in \x are insignificant) the density asymmetry can be maintained at a cost of smaller 
S\x which tends to zero and eventually changes the sign (/Xp > /x„) at Se ~ —1. This marks 
the "turnover" in the gap equation which eventually vanishes for Se > —1. We shall see that 
the regions of the negative values of the deformation parameter are physically not relevant, 
since these states have larger kinetic energies. As in the case of positive Se the effect of finite 
momentum is to suppress the pairing for a large deformation; the combined phase increases 
the pairing only for small perturbations from the asymmetric BCS state in the vicinity of 
the "minimum plateau" {Se ^ 0, P ^ 0). The symmetric part of the quasiparticle spectrum 
(20) controls the location of the maximum (corresponding to Es = 0) and the width of the 
kernel in Eq. (38) in the momentum space through its contribution to the denominator 
of the gap equation. A positive contribution to the average chemical potential effectively 
decreases the density of states (which can be easily seen analytically in the zero temperature 
and weak coupling limits), and therefore, the gap function. 
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FIG. 2. The free energy difference between the superconducting and normal states as a 
function of the Cooper pair momentum P in units of Fermi momentum and the relative deformation 
parameter be. The free energy is normalized to its value in the asymmetric BCS state. The 
parameters are the same as in Fig. 1. 

Inversely, a negative contribution to the average chemical increases the density of states. 
The changes in the Es due to the deformation do not affect qualitatively the results above. 
In summary of Fig. 1, the combination of the LOFF and DFS phases promotes the pairing 
only in the region of small he and P (weak perturbation from the BCS state). For large 
quadrupole deformations the finite momentum is not favorable. Viewed from the LOFF 
phase, for large values of P the quadrupole deformations are disfavored. 

Figure 2 displays the difference between the free energies of the superconducting and 
normal states normalized to its value in the asymmetric BCS state ^JFqo = 5T{P ~ 
0, (5e = 0). The thermal contribution due to the finite-temperature entropy is numerically 
irrelevant for the net budget of free energies of both states. Whether the superconducting 
state is thermodynamically favored depends on the relative magnitude of the potential energy 
of the pair interactions [second term in Eq. (43)] and the difference in the kinetic energies 
of the normal and superconducting states. Since the energy of the pair interactions scales 
as the square of the pairing gap, the shape of the surface closely resembles that of the 
pairing gap in Fig. 1. There are, however, significant quantitative differences due to the 
contribution from the kinetic energy of the quasiparticles. The asymmetric BCS state is the 
stable ground state of the system {J- < 0), however its perturbations for finite 5e and P are 
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unstable towards evolution to lower energy states. For the pure LOFF phase {Se = 0) the 
ground state corresponds to finite momentum P ~ 0.5 (in units of pp). For the pure DFS 
phase (P = 0) there are two minima corresponding to 6e ~ —0.8 and 6e ~ 0.55, i.e., prolate 
and oblate deformations of neutron Fermi spheres, respectively. In general the position of 
the minimum of ST in the Se-P plane (passing through the minima of the limiting cases) 
prefers either large deformations or large finite momenta. The absolute minimum energy 
state corresponds to 6e ~ 0.55 and P = 0; that is, while the LOFF phase is a local minimum 
state, it is generally unstable towards evolution to a pure DFS phase with oblate/prolate 
deformations of neutron/proton Fermi spheres. The effect of the kinetic energy contribution 
is to suppress the large pairing contribution for negative values of the deformation parameter. 
The position of the true minimum of the free energy coincides with the maximum in the 
pairing gap in the positive Se regime. 



B. LOFF phase at finite-temperature 

When the system is in an isospin asymmetric state the effect of the temperature on the 
pairing gap and related characteristics of the superconducting state is twofold: first, at high 
temperatures the pairing is suppressed due to the thermal excitation of the quasiparticle 
states, in analogy to the classical BCS superconductors; second, at low temperatures the 




FIG. 3. The LOFF pairing gap as a function of the Cooper pair momentum and the tempera- 
ture. The parameter values are the same as in Fig. 1. 
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pairing is enhanced by the temperature, since the thermal excitation of the quasiparticle 
states at different Fermi surfaces reduces the blocking effect and increases the phase-space 
overlap. Figure 3 shows the pairing gap in the LOFF state as a function of the temperature 
and finite momentum. For P — the asymmetric BCS state exists between two critical 
temperatures corresponding to the thermal stimulation and suppression of the pairing. For 
constant asymmetry the temperature reduces the antisymmetric part of the quasiparticle 
spectrum [Ea in Eq. (38)], and hence the decoherence between the quasiparticle states on 
different Fermi surfaces. This is responsible for the gap reentrance phenomenon at lower 
critical temperature with increasing temperature. The reduction of the contribution from 
the symmetric part Eg of the spectrum with increasing temperature is responsible for the 
quenching of the superconducting state at the upper critical temperature. At T = 0, but 
for finite P, analogous reentrace of the gap function is seen in Fig. 3. The existence of the 
lower and upper critical momenta can again be understood from the contributions of the 
finite momentum to the quasiparticle spectrum. 




FIG. 4. The free energy difference between the superconducting and normal states SJ^ as a 
function of the Cooper pair momentum P in units of Fermi momentum and temperature. The 
parameters are the same as in Fig. 1. 

In the general case P 7^ and T 7^ the reentrace phenomenon extends in the P — T 
plane: the pairing region of the phase diagram is bounded by the curves for the lower and 
upper critical values Pci{T), Pc2{T), or equivalently Tci{P), Tc2{P). The duality between the 
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effect of finite T on the pairing gap at fixed P to that of finite P at fixed T is evident. In Fig. 
4 we show the difference between the free energies of the superconducting and normal states 
for the LOFF phase. The true ground state here corresponds to the zero-temperature 
hmit for pair momentum P/pp ~ 0.4. The effect of the temperature on the LOFF phase 
depends on the momentum of the pair: for low momenta there is a minimum as a function 
of the temperature (as is also the case for the asymmetric BCS state). For large fixed 
momenta the temperature reduces the pairing gap. The same applies for the perturbations 
from the P = state at fixed temperature: for low temperatures the system prefers finite 
momenta, while for large temperatures the finite momentum only reduces the energy. The 
overall dependence of the free energy on the temperature and the finite momentum reflects 
the dependence of the pairing gap and the potential energy of the pairing interaction on 
these quantities. 

C. DFS phase at finite-temperature 

The features seen in the phase diagram of the DFS phase at finite temperatures are 
analogous to those already discussed for the LOFF phase. Here we reiterate this discussion 
with the emphasis on the qualitative differences of these phases. The pairing gap as a 
function of the temperature and the relative deformation {5e) is shown in Fig. 5. 




FIG. 5. The DFS pairing gap as a function of the relative deformation and the temperature. 
The parameter values are the same as in Fig. 1. 
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The reentracc effect for T = is seen for both signs of Se; the lower critical deformation 
corresponds to the onset of the pairing due to the restoration of the phase-space coherence 
via the deformation of the spherical Fermi surface into quasiellipsoidal form. The upper 
critical Se corresponds to the loss of coherence due to the large shift in the Fermi surfaces 
of the paired particles. The T — Se duahty is seen in the case where deformation vanishes 
(asymmetric BCS state). The lower critical temperature corresponds to the thermal stim- 
ulation of the pair correlations due to the smearing of the Fermi surfaces, while the upper 
critical temperature corresponds to the ordinary thermal suppression of the pairing in the 
BCS-normal state transition. The effect of Se at finite-temperatures can be understood in 
terms of its contribution to the antisymmetric part of the quasiparticle spectrum, in com- 
plete analogy to the LOFF phase. The reentrace effect extends in the plane T — Se, and the 
pairing region is bounded by the critical values Sei(T) and Se2{T), or equivalently Tci{Se) 
and Tc2{Se). Unlike the LOFF phase these curves do not terminate at Se = 0, rather they 
show a continuous crossover from the oblate to the prolate deformation of the Fermi surface 
of the paired particles. For either sign oi Se the temperature thermally stimulates the su- 
perconducting phase for small perturbations from the asymmetric BCS state corresponding 
to the Se — cut; for large deformations the pairing is largest in the zero-temperature limit. 




FIG. 6. The free energy difference between the superconducting and normal states 6J-' as a 
function of the relative deformation and temperature. The parameters are the same as in Fig. 1. 
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The effect of finite Se at fixed temperature is a suppression of tfie pairing in tfie fiigfi- 
temperature limit, wliile at lower temperatures there is a pronoimced maximum in the 
pairing gap as a function of 5e. Finally, the difference between the free energies of the su- 
perconducting and normal states ST for the DFS phase is shown in Fig. 6. There are two 
minima corresponding to the oblate or prolate deformations of the Fermi surfaces of neu- 
trons/protons. The shape of 5 J- is controlled by the interplay between the pairing potential 
energy and the kinetic energy; the contribution from the entropy difference is numerically 
negligible. The minima correspond to large deformations of the order 0.5 which optimize 
the gain in the potential energy over the loss in the kinetic energy due to the deformation; 
the true minimum occurs for positive Se at zero-temperature. While temperature stimu- 
lated pairing correlations are seen for small deviations from the BCS state, this effect is not 
essential, since the minimum of the 5jF surface occurs at large Se. The DFS phase is most 
stable in the zero-temperature limit. 

IV. SUMMARY AND OUTLOOK 

The fermionic condensation in nuclear matter under isospin asymmetry leads to super- 
conducting states which spontaneously break the space symmetries. We studied two realiza- 
tions of such condensates - the LOFF phase, where the Cooper pairs carry finite momentum, 
and the DFS phase where the Fermi surfaces of the paired nuclcons arc deformed, to the 
lowest order into ellipsoidal form. A unifying feature of these phases is the deformation 
of the spectrum of the quasiparticle excitations at the leading (LOFF) and next-to-leading 
(DFS) order in the expansion in Legendre polynomials with respect to the angle between 
the particle momenta and the axis of symmetry breaking. When these deformations are 
treated simultaneously, we find that the pure DFS phase has lower energy than the LOFF 
phase. The combination of these phases still lower the energy of the normal system and 
can exists in a metastable state. In particular, the pure LOFF state could be reahzed in 
the case where the unpaired state features some intrinsic spatial asymmetries (e.g., due to 
the lattice in the metallic superconductors). Of course, the axis of the symmetry breaking 
for the LOFF and DFS phases need not be the same and, more generally, the patterns of 
symmetry breaking could be more complicated than those studied here (e.g., higher order 
multipole deformations of the Fermi surfaces) . 

The phase diagrams of the LOFF and DFS phases in the temperature/deformation plane 
show a number of unusual features: (i) two critical temperatures for the onset of the su- 
perconducting state (the reentrace phenomenon). The upper critical temperature is the 
analog of the BCS temperature for the superconducting-normal phase transition. The lower 
critical temperature is specific to asymmetric systems and results from the thermally stim- 
ulated phase-space overlap between the pairing states, (ii) There exists duality between 
the temperature effects and the finite-momentum or deformation effects. For example, at 
low temperatures the LOFF phase is characterized by two critical momenta of the Cooper 
pairs at which the superconducting phase sets in. (iii) The free energy of the combined 
LOFF and DFS phases shows local minima corresponding to negative and positive relative 
deformations at zero momentum of the pairs. The lowest energy corresponds to a phase 
with oblate neutron and prolate proton Fermi surfaces {Se is positive). 

The realization of the ^5*1 — ^ Di pairing in bulk nuclear matter depends on a number 
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of factors: (i) the ratio of the pairing gap to the average chemical potential, which controls 
the dcnsity-temperature-dcpendent critical asymmetries at which the pairing vanishes; this 
ratio depends sensitively on the quasiparticle renormalization effects which need further 
study in this context, (ii) A second factor involves the asymmetries imphed by the /? 
equilibrium of the matter with leptons in the stellar matter; while the asymmetries derived 
for the normal state may serve as order of magnitude estimates, the LOFF and DFS phases 
under f] equilibrium might require asymmetries largely different from those deduced for 
the unpaired system, (iii) The matter could prefer a state where the superconducting and 
normal phases coexist in a mixed phase, with superconducting "islands" immersed in the 
unpaired matter; the analogy to the type-I superconductors in the magnetic fields is a useful 
counterpart of this type of ordering, (iv) The extreme low-density asymmetric nuclear 
matter can support pairing in the ^5*1 — ^ Di channel at much larger asymmetries than the 
matter, say, at the nuclear saturation density; however, in this strong coupling regime one 
deals with a Bose-Einstein condensate of deuterons coexisting with a low-density neutron 
gas [25] and one needs an understanding of the crossover from the LOFF and DFS phases 
to their Bose-Einstein condensed counterparts. 

We now briefly comment on possible realization of the phases with broken space sym- 
metries in other nuclear systems. The neutron-proton pairing in finite nuclei [48-58] is an 
obvious candidate. The odd-number projected BCS equations contain the decoherence ef- 
fects in the gap equation introduced by the extra odd particle. The resulting blocking effect 
is quite analogous to the phase decoherence due to the imposed neutron-proton asymmetry 
in the bulk nuclear matter. The problem of realization of the LOFF phase is facilitated by 
the fact that the pairing instability occurs for total momenta P ^ pp ^ ^/R, where R is the 
radius of the nucleus. Therefore, for already medium size nuclei the influence of the surface 
on the condensation energy should not be large and the solutions found for the interior of 
the nucleus should go over to those corresponding to infinite systems. The translational 
symmetry is not broken by the DFS phase, therefore the finite size effects (while being gen- 
erally of the same order of magnitude as for the LOFF phase) would affect the structure 
of the condensate to the same extent as those for the case of the pairing with undeformed 
Fermi surfaces. 

Quite generally, the phases with broken space symmetries arise in the physical situations 
where the time-reversal symmetry between the paired fermion states is broken. This is the 
case, for example, in finite nuclei under rotation. The ^5*0 pairing among neutrons and/or 
protons in such nuclei can support the LOFF or DFS like states. Similarly, the magnetic field 
breaks the time-reversal invariance due to the Pauli paramagnetism of neutrons and protons. 
The LOFF and DFS phases can arise for the case of the ^Sq pairing in strong magnetic fields, 
for example, in neutron star crusts. These phases will emerge when the interaction energy of 
a neutron spin with the magnetic field fisB, where /jb is the anomalous magnetic moment of 
a neutron, is of the order of the pairing gap in the zero field limit. This condition is satisfied 
for the field strength of the order 10^^-10^^ G (depending on the value of the pairing gap) 
which could occur in the highly magnetized subpopulation of neutron stars (magnetars). 

As is well know, the breaking of the continuous symmetries in the infinite systems leads to 
collective excitations with vanishing minimal frequency (Goldstone's theorem). The LOFF 
phase breaks the translational and rotational symmetries; the existence of the preferred 
direction in the DFS phase breaks the rotational symmetry of the system. The breaking 
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of continuous symmetries in the ground of these phases, therefore, imphes new collective 
bosonic modes in asymmetric nuclear matter. In finite systems the Goldstone theorem 
implies collective modes which tend to zero as a certain power of the radius of the system. 
The observation of the Goldstone modes associated with the LOFF and/or DFS phases could 
provide information on the true ground state structure of the condensate under laboratory 
conditions, for example, in finite nuclei. 
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